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1. Introduction 

One of the most important quantities in the low energy kaon physics is the form factor / + at 
zero momentum transfer, which provides a theoretical input for the determination of the Kobayashi- 
Maskawa matrix element \V US \ through the decay. 

For this form factor it is known that 

1. The leading contribution is unity, which is exact in the SU(3) limit, 

2. There is no contribution of 0{m s — m) due to the Ademollo-Gatto theorem where m = (m u + 
m d )/2. 

In 1984, Leutwyler and Roos estimated this form factor analytically [|T]] and their result /+(0) = 
0.961 ±0.008 has been used as the standard in the phenomenological analysis. Though the leading 
correction was determined unambiguously using chiral perturbation theory, the next-to-leading 
order correction was estimated using a model of the wave function of the pseudoscalar meson, with 
which only a crude error estimate is available. The numerical simulation of lattice QCD could 
help this situation. Indeed, a quenched result appears recently [Q] and several groups carry out the 
unquenched calculation with different lattice setups [||]. 

Our calculation is performed on two-flavor QCD gauge configurations generated with the non- 
perturbatively 0(a)-improved Wilson fermion by the JLQCD Collaboration The lattice size is 
20 3 x 48 at j8 = 5.2 and 12,000 HMC trajectories are accumulated. Five sea quark masses corre- 
sponding to mps/my ~ 0.8 — 0.6 enable us to make a careful study of the quark mass dependence 
of the form factors. Especially, we focus on the consistency between the lattice data and the ChPT 
predictions. The pion form factors are obtained as a by-product of the calculation. The results are 
presented in the poster session by S. Hashimoto 



2. Kaon form factors 

The kaon semileptonic decay form factors f+{q 2 ) and f-{q 2 ) are defined as 

(K(p)\Vv\K(k))=U(q 2 )(k + p) tl +f_(q 2 )(k-p)t 1 , (2.1) 

where q 2 = (k — p) 2 and is the vector part of the weak current. Reliable theoretical estimate of 
/+(0) is necessary for the precise determination of the Kobayashi-Maskawa matrix element \V US \ 
from the decay. Our goal is to calculate the form factor in a few percent accuracy from the first 
principle of QCD. To do that, careful study of SU(3) breaking effects and the chiral extrapolation 
are vital. 

The scalar form factor fo(q 2 ) can be expressed as a linear combination of f+{q 2 ) and /_ (q 2 ) 

Mq 2 )=f + (q 2 )+ / 2 f-(q 2 ), (2-2) 

m K~ m K 

and the ratio of /_ (q 2 ) and f+(q 2 ) is denoted as t, (q 2 ) 
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Figure 1: Double ratio R K ^ K (t) for the scalar form factor /o(^„ ax ) at K sea = 0.1355 



3. Lattice calculation 



Our calculation proceeds in three steps, and the form factor / + (0) is expressed by a product of 
three factors 



/+(0) 



f+{0) = fo(q 2 max ) X 



1 + 5(0) 



m K -m n 
niK + mjz 



1 



Mllua) 



1 + 5(0) 



m K -m n 
m K + m„ 



(3.1) 



In the first step, the scalar form factor at the maximum momentum transfer squared q 2 max = {rtiK — 
m n ) 2 is estimated and the interpolation to q 2 = is done in the second step. The unnecessary 
contribution from £ (0) is subtracted in the last step. A double ratio of the three-point functions is 
used to obtain each factor. We explain the details of each step in the following subsections. 

3.1 Double ratio I 

The scalar form factor at the maximum momentum transfer squared /oGziLa) * s obtained from 
the double ratio of three-point functions 



R 



(0 



C %ViK (t,T/2;d,d)C mn (t,T/2;d,d) 



Cx V4n {t,T/2;0,0)C mK (t,T/2;0,0) 
where the three-point function is given by 



m K + m % , 2 



M 



, o«a<r/2, (3.2) 



(3.3) 



and O n , Ok is the interpolating operator of pion and kaon. This double ratio is essentially the same 
as that used before in the calculation of the B — » Dlv form factor by the Fermilab group [0]. 
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Figure 2: The ratio of form factor /+(?)//+ ( f o) (left) and £, (t) (right) as a function of momentum transfer 
squared t = q 2 at K sea = 0.1355. A notation f + (t) = f+(t) [l + % (O ^gjg 



is used here. 



A typical time dependence of the double ratio R K ^ n (t) at the lightest sea quark mass K sea = 
0.1355 is plotted in Figure [[]. Smeared kaon andpion sources are set at time t = 0, T/2, respectively. 
The vector current is inserted at t. Because not only the bulk of the statistical fluctuation but 
the multiplicative renormalization factors and other systematic errors cancel, we can determine 
fo(qliax) w i tn an accuracy better than one percent. 

3.2 Double ratio II 

The second double ratio with the pion two-point function C %n {t; p) 
C7tv 4 K(t,T/2;p,p) 



C mK (tJ/2;0,0) 
Cmz{t\p) 



mK + E n (p) 


r ^-iwffli 

m K +E K (p) 


m K + m n , 2 

J+Klmax) 


mK + m n 



< t <C T/2, (3.4) 



is used to interpolate the form factor to q 2 = 0. This is the second factor in ( |3.1[ ), because fo(q 

f+ (^max ) 1 + <■> (?Lr ) 



Imax/ 



max) m K +m n 



. In this case, the momentum is inserted at t and T/2 so that the pion 
has a finite momentum and the kaon is at rest. The double ratio ( |3.4[ ) is statistically noisier than 
the double ratio ([T2|) in the previous subsection. The data corresponding to the pion momenta 
p = 2n/20 x (±1,0,0), (±1, ±1,0) and (±1,±1,±1) are obtained and the quadratic function is 
used to interpolate to q 2 = 0, which is shown in Figure || (left). 

3.3 Double ratio III 

The last double ratio is for the determination of ^ (0) 

C nViK {t, T/2;p,p)C mjl (t, T/2; p,p) 1 - § (q 2 ) 



C nV4K (t,T /2;p,p)C KVi% {t,T /2;p,p) m K + E n {p) | ^^ m K -E n {p) 



0<?«T/2, (3.5) 



m n +E n {p) m n +E n (p) 

where we have to measure the spatial component of the matrix element. The extrapolation to q 2 = 
is done by assuming the linear dependence on q 2 , which is shown in Figure ^| (right). We see that 
the q 2 dependence of the E, is very weak and independent of the strange quark mass. 
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Figure 3: Chiral extrapolation of the form factor /+ (0) with quadratic function (left) and the one-loop ChPT 
formula plus a quadratic function (right). 



4. Result for / + (0) 

Following the method explained in the previous subsections, the form factor /+ (0) is estimated 
as a function of the up/down and strange quark masses. In order to obtain the form factor at the 
physical pion and kaon masses, we first use a simple fitting function 



/+(0) = 1 - (c +c\ [(r Q m K ) 2 + (r m„) ])[(r m K ) 



(ro«% 



\2-i2 



(4.1) 



where the meson masses are normalized by the Sommer scale r$. The result is shown in Figure ||. 
Our data are well represented by function Q4.1[ ) and the preliminary result /+(0) = 0.967(6) is con- 
sistent with Leutwyler-Roos's value 0.961(8) [JJ], and the quenched result 0.960(9) by Becirevic 
et al. [@] 

For small pion and kaon masses the ChPT predicts the mass dependence of the form factor 



2raprag 



In- 



nr, 



/+(<>) = l + |/fe*(0) + |/fr„(0), *fi(°) = -i28^ 
We try to fit the data with 

/+(0) = 1 + 3 -H KK {0) + 3 -H Kn (0) - (c + cy [(r m K ) 2 + (r m K f])[(r m K f - (r Q m K f] 



2 2 2 

nip — niQ mp 



2i2 



(4.2) 



(4.3) 



At this order the one-loop ChPT formula has no tunable parameters. The chiral logarithm is signif- 
icant only in the region where m\ <C m 2 K , while the data in the region 1/2 < m 2 n /m 2 K < 2 are well 



approximated by the quadratic form fl4.ip. The extrapolated value at physical pion and kaon mass 
is affected by the chiral logarithm. The result is /+(0) = 0.952(6). 



5. Charge radius 

The charge radius (r 2 ) parametrizes the slope of the form factor near q 2 = 

1 



(5.1) 
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The parameter A = {r 2 )m 2 n /6 is often used in the literature. 

For the vector form factor, our result extrapolated with the one-loop ChPT plus a quadratic 
function shown in Figure |] is (r 2 )y K = 0.26(3) fm 2 . It corresponds to A+ = 0.021 (2), which is sig- 
nificantly smaller than the experimental value 0.0278(7) or the quenched lattice result 0.026(2) 
by Becirevic et al. [Q] Our result for the charge radius of the scalar form factor (Figure |]) is 
(r 2 )^ 71 = 0.37(6) fm 2 , which corresponds to Ao = 0.031(5) that overshoots the experimental value 
0.0174(22) or the quenched lattice result 0.012(2). 

Our calculation shows that the kaon form factors can be calculated on the lattice with good 
precision using the double ratio method. Further study on systematic errors, especially scaling 
violation effects, is needed to obtain the definitive value. 
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